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OO ■ Abstract 

Qs, '.- leading order (NLO) QCD calculations. The dependence of the theoretical predictions on the structure 
On ■ functions, and on the theoretical scales is investigated. It is shown that the data cannot be simulta- 
(~| • neously fitted with a single set of structure functions and scales. However, it is argued that as long as 
Qh, one restricts the data/theory comparison to the xt range where the theory is reliable, i.e. stable with 
O il respect to the scale variation, there is no need to introduce an additional primordial /c^ dependence 
^ ■ except for Be data. Finally a precise determination of the strong coupling constant, q^, is performed 
• • , from the direct photon production cross sections obtained in high statistics pp and pp collisions at the 
.^ CERN SPS (UA6) by a NLO QCD analysis. 
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1 Critical study of photon production in pp and pBe interactions 

Inclusive photon production including the direct and the bremsstrahlung contributions, calculated at 
NLO in QCD [|l|, is compared to experimental data in the range v^=23 GeV to 63 GeV as a function 
of XT = 2pT/\fs. The following experimental data are used: the pp data from WA70 and from 
ISR R806 §, RllO @, AFS [|] and R108 § as weh as the new fixed target data from E706 § {pBe) 
and UA6 ||8[ {pp and pp). The main results of this study are summarized below. They concern the 
theoretical prediction uncertainty and possible incompatibilities between data sets. 

The perturbative NLO calculations involve three arbitrary scales: the factorization scale M, the 
renormalization scale /U and the fragmentation scale Mp. As those scales are unphysical, the theory 
can be considered reliable only in the region of the phase space where the predictions are stable with 
respect to the scale variations. For a large range of Mp values, the theory stability versus M and /i 
is reasonable for Xj, > .35 at ^/s=23 GeV, Xj, > .32 at y^=31.6 GeV, Xj. > .26 at y^=38.8 GeV 
and Xj, > .16 at ^/s=63 GeV. The theoretical uncertainty due to the scales is approximated by the 
difference in the predictions between M = fi = Mp = Pt/'^ and M = ^ = Mp = Pt/2. 

The ratio data/theory is shown in Fig. [l| with all scales set to pt/2 and using the parton distribution 
functions CTEQ4M |^]. A large discrepancy can be seen between the E706 and the other data sets. 
Although there is an approximate agreement with theory for the WA70, UA6 and ISR data sets, the 
E706 data are underestimated by the theory by a factor 2 in the region of theory stability and up to 
4 at low XT- In Fig. g the same experimental data are compared with the theory with all scales set to 
pr/3 and using the parton distribution functions MRS98-2 [|lO|. Here too the theory underestimates 
the E706 data by 50% in the region of theory stability and by up to 100% at low xt- 

Clearly, the full set of available data cannot be fitted by the theory whatever the scales or the struc- 
ture functions. Restricting the comparison to the xt range of theory stability reduces the discrepancy 
to a mere global scaling of the E706 data. 

It has been proposed to introduce a parton-fcy broadening effect 0, |l^] to account for the low 
XT dependence and for the scaling effect seen by E706. This effect would destroy the theoretical 
agreement of all other experiments particularly in xt region of theory stability. 



A NLO QCD calculation has been performed for vr'^ production [12|. Setting all scales to pt/2, 
the low XT rise observed in the direct photon case (Fig. |l]) has almost disappeared (Fig.^. Moreover, 
when all scales are set to pt/3, a common rescaling of the theory by a factor 1.5 it 0.2 is enough to 
accommodate all fixed target data. This normalization may not be unrealistic as the absolute normal- 
ization of the theory is not precise, due for example to present uncertainties in the vr fragmentation. 
The fixed target 7r° data appear in better agreement then the photon data. 



2 Determination of a. from UA6 data [[13 



The measurement in the same experiment of the production of direct photon in pp and pp collisions 
[^ allows a clean isolation of the annihilation process {qq — > 73) from the difference a{pp — > 7X) 
—a{pp -^ jX). Knowing the quark distributions, this process provides a direct handle on AiiL and 
therefore ag. Quark distributions vs A^ have been obtained from fits to BCDMS |[l^ deep inelastic 
scattering data. Therefore A^ remains the only free parameter in the fit of the UA6 data to the 
NLO theoretical calculations. The results are shown in Fig. ^ for various choices of scales. 



For optimized scales |jT5|, the best x^ gives A^ of 210 ± 22 (stat.) ± 44 (syst.) MeV. We conser- 
vatively assign an overall theoretical uncertainty of ^35 MeV, dominated by uncertainty in the fi and 
M scale selection. 

This leads to as(M|) = 0.1112 ± 0.0016 (stat.) ± 0.0033 (syst.) laooH (theo.) in good agreement 
with the world average value [^] and with error bars comparable to those achieved in deep inelastic 
scattering |17, O]. Recent improved NLO computation for pp including resummation shows small 



corrections for small scales as used here and furthermore reduces the scale dependence as expected 
[p^]. Hence these new calculations should not change significantly the Os value. 
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Figure 1: Ratio data/theory for direct photon production using CTEQ4M distribution functions. All 
scales are set to pr/2- The arrows note the x^ ranges where the perturbative predictions are reasonably 
stable vs variation of the scales n and M. Statistical errors are shown as full lines, statistical and 
systematic errors added in quadrature are shown as dashed lines. 
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Figure 2: Ratio data/theory for direct photon production using MRS98-2 distribution functions. All 
scales are pt/3. 
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Figure 3: Ratio data/theory for vr*^ production using CTEQ4M parton distribution functions. All 
scales are set to pt/2- The E706 pBe data with beams of 530 GeV/c and 800 GeV/c as well as the 
E706 7r~Be data with a beam of 515 GeV/c (not discussed here) are shown. 
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Figure 4: a) The x^ between the theoretical predictions and the measured cross section difference 
g-ipp^7 ) g-(pp^7 J sunimed over ten pri bins as a function of A^ for various choices of scales. 
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b) Best value of ALL as a function of the parameter c defining the scales fjr = M = Mp = cpj^. The 
error bars are statistical only. 



